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Abstract 

After new publications of HI data for the diffractive photoproduction of dijets, which overlap 
with the earlier published HI data and the recently published data of the ZEUS collaboration, 
have appeared, we have recalculated the cross sections for this process in next-to-leading order 
(NLO) of perturbative QCD to see whether they can be interpreted consistently. The results of 

> ■ 

these calculations are compared to the data of both collaborations. We find that the NLO cross 



sections disagree with the data, showing that factorization breaking occurs at that order. If direct 



and resolved contributions are both suppressed by the same amount, the global suppression factor 

o 

depends on the transverse-energy cut. However, by suppressing only the resolved contribution, 
also reasonably good agreement with all the data is found with a suppression factor independent 

•rH . 

^ ' of the transverse-energy cut. 
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FIG. 1: Diffractive scattering process ep — > eXY, where the hadronic systems X and Y are 
separated by the largest rapidity gap in the final state. 

I. INTRODUCTION 



At high-energy colliders such as the ep collider HERA at DESY and the pp collider 
Tevatron at Fermilab, diffractive processes are known to constitute an important fraction of 
all scattering events. These events are defined experimentally by the presence of a forward- 
going hadronic system Y with four- momentum py, low mass My (typically a proton that 
remained intact or a low-lying nucleon resonance), small four- momentum transfer t = (p — 
Py) 2 , and small longitudinal momentum transfer xjp = q(p — Py)/{(1P) from the incoming 
proton with four-momentum p to the central hadronic system X (see Fig. [T] for the case of 
ep — > eXY). Experimentally a large rapidity gap separates the hadronic system X with 
invariant mass Mx from the final-state system Y. 

Theoretically diffractive interactions are described in the framework of Regge theory [l| as 
the exchange of a trajectory with vacuum quantum numbers, the Pomeron (IP) trajectory. 
Then the object exchanged between the systems X and Y, as indicated in Fig. (TJ is the 
Pomeron (or additional lower-lying Regge poles), and the upper vertex of ep — > eXY, i.e. 



2 



eP — > eX, can be interpreted as deep-inelastic scattering (DIS) on the Pomeron target 
for the case that the virtuality of the exchanged photon Q 2 = —q 2 is sufficiently large. In 
analogy to DIS on a proton target, ep — > eX, the cross section for the process eP — > eX in 
the DIS region can be expressed as the convolution of partonic cross sections and universal 
parton distribution functions (PDFs) of the Pomeron. The partonic cross sections are the 
same as for ep DIS. The Pomeron PDFs are usually multiplied with vertex functions for 
the lower vertex in Fig. [H yielding the diffractive parton distribution functions (DPDFs). 
The Q 2 evolution of the DPDFs is calculated with the usual DGLAP [2] evolution equations 
known from ep — > eX DIS. Except for the Q 2 evolution, the DPDFs can not be calculated 
in the framework of perturbative QCD and must be determined from experiment. Such 
DPDFs 3j-[6| have been obtained from the HERA inclusive measurements of the diffractive 
structure function , defined in analogy with the proton structure function F^. 

The presence of a hard scale such as the squared photon virtuality Q 2 = —q 2 in deep- 
inelastic scattering or a large transverse jet energy E^ 1 in the photon-proton center-of- 
momentum frame should then allow for calculations of partonic cross sections for the central 
system X using perturbative QCD. Such diffractive processes with the presence of a hard 
scale are usually called hard diffractive processes. The central problem in hard diffraction 
is the problem of QCD factorization, i.e. the question whether diffractive cross sections are 
factorisable into universal diffractive parton density functions and partonic cross sections. 

For DIS processes, factorization has indeed been proven to hold [7| , and DPDFs have been 
extracted at low and intermediate Q 2 (si, 4| from high-precision inclusive measurements 
of the process ep — > eXY using the usual DGLAP evolution equations. The proof of 
the factorization formula, usually referred to as the validity of QCD factorization in hard 
diffraction, also appears to be valid for the direct part of photoproduction Q 2 ~ or low- 
Q 2 electroproduction of jets Q]. However, factorization does not hold for hard processes 
in diffractive hadron-hadron scattering. The problem is that soft interactions between the 
ingoing hadrons and their remnants occur in both the initial and final states. This agrees 
with experimental measurements at the Tevatron [8|. Predictions of diffractive dijet cross 



sections for collisions as measured by CDF using DPDFs determined a few years ago 
and more recently [41 by the HI collaboration at HERA overestimate the measured cross 

n □ 

section by up to an order of mag nitude [8|,H0[. This suppression of the CDF cross section 
can be explained by the rescattering of the two incoming hadron beams which, by creating 
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FIG. 2: Diffractive production of dijets with invariant mass Myi in direct (left) and resolved (right) 
photon-pomeron collisions, leading to the production of one or two additional remnant jets. 



additional hadrons, destroy the rapidity gap ll) . Jet production with real photons involves 
direct interactions of the photon with quarks or gluons from the proton (or in our case from 
the pomeron) as well as resolved photon contributions, leading to parton-parton interactions 
and an additional remnant jet coming from the photon as reviewed in 12] (see Fig. |2]). For 
the direct interactions, factorization is expected to be valid as in the case of DIS, whereas 
we expect it to fail for the resolved process as in hadron-hadron scattering. For this part of 
photoproduction one would therefore expect a similar suppression factor due to rescattering 
effects of the ingoing partons. Introducing vector-meson dominance photon fluctuations, 
such a suppression by about a factor of three for resolved photoproduction at HERA was 
predicted 13]. 

On the experimental side, the first measurements of dijet cross sections in diffractive 
photoproduction have been presented by the HI collaboration as contributions to two con- 
ferences 14] . The kinematic range for these data were Q 2 < 0.01 GeV 2 , xjp < 0.03, Erf 1 > 5 
GeV, E 3 t > 4 GeV and 165 < W < 240 GeV, where jets were identified using the inclu- 
sive A/r-cluster algorithm. The measured cross sections as a function of x° bs and z^ s were 
compared to leading-order (LO) QCD predictions, using the RAPGAP Monte Carlo model 
0. For the DPDFs the LO 'HI 2002 fit' was used fl. It was found that these two cross 
sections were well described by the predictions in normalization and shape over the whole 
range of x° bs and z°p , showing no breakdown of factorization in either the resolved or in the 
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direct photoproduction. In addition, normalized cross sections as a function of various other 
variables were compared to the predictions with the result that all measured distributions 
were in good agreement. 

Subsequently we calculated the next-to-leading order (NLO) corrections for the cross 
section of diffractive dijet production using the same kinematic cuts and with the same 
DPDFs as in the first HI analysis 14|on the basis of our previous work on NLO corrections 
for inclusive direct [16] and resolved 17J dijet photoproduction. While at LO good agreement 
with the HI data [l4| was found, consistent with the finding in the HI analysis 14J, it was 



found that the NLO corrections increase the cross section significantly 
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19| and require a 



suppression factor of the order of R = 0.5. Since on theoretical grounds only a suppression 
of the resolved cross section would be acceptable, we demonstrated in Q that by 
multiplying the resolved cross section with a suppression factor of R = 0.34, reasonably 
good agreement with the preliminary HI data 14J could be achieved. This value for the 
suppression factor turned out to be in good agreement with the prediction of {"Lsl ]. 

The first experimental data from the ZEUS collaboration were presented at the DIS work- 
shop in 2004 [20(1 . The dijet cross sections were obtained in the kinematic range Q 2 < 1 
GeV 2 , xjp < 0.025 and E^* > 7.5 (6.5) GeV. For these kinematic constraints NLO calcu- 
lations were not available in 2004. So, the measurements were compared to LO calculations, 



unfortunately with previous HI DPDFs 



21] with the result that good agreement in the shape 



was achieved. However, the normalization was off by a factor of 0.6, which was attributed 
later to the older DPDF input 22 ]. so that the HI and ZEUS results were consistent with 
each other. The situation concerning the agreement of HI and ZEUS data and the influence 
of NLO corrections improved already considerably in the fall of 2004. At the ICHEP 2004 
both collaborations presented their new data and compared them with NLO cross section 
calculations. HI compared their data with the predictions from the program of Frixione 
23 ] and ZEUS with our calculations along the lines of Q, [3], where now only the different 
kinematic cuts of the ZEUS analysis had to be incorporated. Both collaborations also used 
the same DPDFs, namely the 'HI 2002 fit' 9]. The conclusion from the comparison of the 
respective NLO calculations with the HI [24J and ZEUS [25] measurements were very sim- 
ilar. Both collaborations observed from their data that good agreement was achieved with 
the global suppression (of direct and resolved contribution) by a factor 0.5. Concerning the 
model with suppression of the resolved contribution only, the HI collaboration concluded 
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that this model was disfavored compared to the global suppression model. This result was 
obtained by comparing the differential cross sections da/dx° bs and da/dy to a calculation, in 
which the partonic cross sections were suppressed for x° bs < 0.9 by a factor of R = 0.34. In 
the ZEUS contribution to the ICHEP 2004, differential cross sections of several observables 
had been shown and compared with our calculations. All five comparisons, namely for y, 
xjp, Zp s , E 3 ^ tl and rf etl , showed very good agreement both in shape and normalization of 
the cross section with the resolved photon suppression (R = 0.34). The only exception was 
the comparison of da/dx° bs , which did not show agreement with the resolved photon sup- 
pression only. Here the suppression factor varied as a function of x° bs between 0.5 and 0.75 
and the description of the shape was better without any suppression, which would signal 
a global suppression. We emphasize that at this stage of the analysis at the end of 2004, 
both the HI |24j and ZEUS 25] collaborations showed that their data were consistent with 
a global suppression of about a factor of two against the sum of direct and resolved NLO 
QCD predictions. In addition, in the HI contribution 24] it was claimed explicitly that 
there was evidence for factorization breaking also in direct photoproduction. In 2005 the 
ZEUS collaboration presented a more detailed comparison by dividing their data sample into 
two subsets, x° bs > 0.75 and x° bs < 0.75, in order to be more sensitive to the suppression 
of the direct (x° bs > 0.75) and the resolved (x° bs < 0.75) components. Together with the 
result of HI the overall conclusion at the DIS 2005 workshop was the same as at the end of 
2004, namely, that good agreement was achieved with the global suppression of 0.5, while a 
suppression of only the resolved contribution at NLO was disfavored by the data |26l |. 

The analysis of the ZEUS data with respect to the samples enriched in direct and resolved 
processes was continued in a contribution to the Uppsala Lepton-Photon conference in 2005 



23| . For x obs > 0.75 the NLO predictions gave a good description of the shape of the 



measured cross section, although the absolute normalization was a factor of two above the 
data. For x obs < 0.75 the NLO calculations were again above the data when no suppression 
(R = 1) was applied and below the data by a factor of two when a suppression with 
R = 0.34 was applied to the resolved photon processes. The ratio of the resolved-enriched 
to the direct-enriched samples was reasonably well reproduced by the NLO predictions with 
R = 1, indicating that a suppression of the resolved sample with respect to the direct sample 
was not seen in any particular kinematic region. This agreed with the earlier findings that 
a uniform suppression for both resolved and direct process gives a better description of the 
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data. Of course, all these conclusions relied on the fact, that the DPDFs as evaluated by 
HI 9] are really the correct ones. The analysis in {2?! was based on the largest selection of 
variables so far, namely y, xjp, Mx, z^ s , and rf etl . 

The conclusions above concerning the global overall suppression versus a suppression in 
the resolved contribution only based solely on preliminary data from HI and ZEUS and 
the preliminary 'HI 2002 fit' jsj for the DPDFs remained also after 2005 until the final 
publications of the HI and ZEUS analysis appeared in 2007. The comparison between the 
final experimental results and the NLO theory used the new and final DPDFs constructed by 
the HI collaboration [4]. This analysis was based on the larger sample of the years 1997-2000 
as compared to the previous published PDF sets. In [4J two NLO fits, 'HI 2006 fit A' and 'HI 
2006 fit B' were presented, which both give a good description of inclusive diffraction. These 
two sets of PDFs differ mainly in the gluon density at large fractional parton momentum, 
which is poorly constrained by the inclusive diffractive scattering data, since there is no direct 
coupling of the photon to gluons, so that the gluon density is constrained only through the 
evolution. The gluon density of fit A is peaked at the starting scale at large fractional 
momentum, whereas the fit B is flat in that region. 



The differential cross sections as measured in diffractive photoproduction by HI 28] were 
compared with the NLO predictions obtained with the Frixione program 23), interfaced to 
the 'HI 2006 fit B' DPDFs. In this publication [28], the conclusions deduced earlier from 
the comparison with the preliminary data and the preliminary 'HI 2002 fit' j^] are fully 
confirmed, now also with the new DPDF fits [4j. In particular, the global suppression is 
obtained, independent of the DPDF fits used, i.e. fit A or fit B, by considering the ratio of 
measured dijet cross section to NLO predictions in photoproduction in relation to the same 
ratio in DIS. In this comparison the value of the suppression is 0.5 ± 0.1. In addition, by 
using the overall suppression factor 0.5, HI obtained a good description of all the measured 
distributions in the variables z^ s , x° bs , xjp, W, E^ a , ff°^ t , |A^ et | and M 12 interfaced with 
the 'HI 2006 fit B' DPDFs and taking into account hadronization corrections 28J]. Finally, 
the HI collaboration investigated how well the data are describable under the assumption 
that in the NLO calculation the cross section for x° bs > 0.9 is not suppressed. The best 
agreement in a fit was obtained for a suppression factor 0.44 for the NLO calculation with 
x° bs < 0.9, based on fitting the distributions for x° bs , W, ff-f t and EjT a . In this comparison 
they found disagreement for the largest x° bs -bin and the lowest f] lab t -hm (which are related), 
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but better agreement in the i?^ e<1 -distribution. In 28j this leads to the statement, that the 
assumption that the direct cross section obeys factorization is strongly disfavored by their 
analysis. In total, it is obvious that in the final HI analysis {28 1 a global suppression in 
diffractive dijet photoproduction is clearly established. 

Just recently also the ZEUS collaboration presented their final result on diffractive dijet 
photoproduction j^]. As in their preliminary analysis, the two jets with the highest trans- 
verse energies E 3 ^ 1 were required to satisfy E^ n ^ > 7.5 (6.5) GeV, which is higher than in 



the HI analysis with E^ tl ^ > 5 (4) GeV [28|. ZEUS compared their measurements with 



the NLO predictions for diffractive photoproduction of dijets based on our program [19]. 
Three sets of DPDFs were used, the ZEUS LPS fit, determined from a NLO analysis of 
inclusive diffraction and diffractive charm-production data js], and the two HI fits, 'HI 2006 



fits A,B' |4|. The NLO results obtained with the two HI fits were scaled down by a factor 
of 0.87 |4J since the HI measurements used to derive the DPDFs include low-mass proton 
dissociative processes with My < 1.6 GeV, which increases the photon-diffractive cross sec- 
tion by 1.15i ; 08 as compared to the pure proton final state as corrected to in the ZEUS 
analysis. The comparison of the measured cross sections and the theoretical predictions was 
based on the differential cross sections in the variables y, My, x P , zp, E J T et \ r]i etl and 
x° bs . The data were reasonably well described in their shape as a function of these variables 
and lay systematically below the predictions. The predictions for the three DPDFs differed 
appreciably. The cross sections for the 'HI 2006 fit A' ('HI 2006 fit B') were the highest 
(lowest), and the one for the ZEUS LPS fit lay between the two, but nearer to the fit A than 
the fit B predictions. For da/dx° bs ZEUS also showed the ratio of the data and the NLO 
predictions using the ZEUS LPS fit. The ratio was consistent with a suppression factor of 0.7 
independent of x° bs . This suppression factor depended on the DPDFs and ranged between 
0.6 ('HI 2006 fit A') and 0.9 ('HI 2006 fit B'). Taking into account the scale dependence of 
the theoretical predictions the ratio was outside the theoretical uncertainty for the ZEUS 
LPS fit and the 'HI 2006 fit A', but not for the 'HI 2006 fit B'. In their conclusions the 
authors of the ZEUS analysis [29] made the statement that the NLO calculations tend to 
overestimate the measured cross section, which would mean that a suppression is present. 
Unfortunately, however, they continued, that, within the large uncertainties of the NLO 
calculations, the data were compatible with the QCD calculations, i.e. with no suppression. 
Such a statement clearly contradicts the result of the HI collaboration 28( and casts 
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doubts on the correctness of the HI analysis. The authors of [29j attribute this discrepancy 



to the fact that the HI measurements 1 28] were carried out in a lower and a higher xw 
range than those in the ZEUS study 29j. Besides the different E^ 1 and xp regions in 28 1 
and 



29] the two measurements suffer also from different experimental cuts of some other 
variables which makes it difficult to compare the two data sets directly (note also the lower 



center-of-mass energ 



second analysis 



30 



V 



or the HI data). With this in mind the HI collaboration has done a 



31], in which most of the experimental cuts are taken as in the ZEUS 



29] analysis, i.e. the cuts on xp, r]^ et1 ^ and on E 3 ^ tl<y2 \ In addition they analyzed data 



sets also with the lower E^ 1 cut, namely E^ n ^ > 5 (4) GeV and with xjp < 0.03 as in the 



previous HI dijet analysis 28]. Starting from these recent data 30|, |3l| we have performed 



Uj for the new HI [30|, |31| and 



a new calculation of the NLO cross sections on the basis of 
the latest ZEUS 29J analyses with the same DPDFs as input, in order to see whether we 



and the ZEUS 



29] 



can confirm the different conclusions obtained from the older HI 
measurements. In this new comparison between the experimental and the theoretical results 
we shall concentrate on using the 'HI 2006 fit B' as DPDF input, since it leads to smaller 
NLO cross sections than the DPDFs based on the 'HI 2006 fit A' or the ZEUS LPS fit. 

In section 2 we shall present, after defining the complete list of cuts on the experimental 
variables and giving all the input used in the cross section calculations, the comparison 



with the new HI experimental data 



30 



3l| . In this comparison we shall concentrate on 



the main question, whether there is a suppression in the photoproduction data at all. In 
addition we shall investigate also whether a reasonable description of the data is possible 
with suppression of the resolved cross section only, as we studied it already in our previous 
work in 2004 



18 



19] . In section 3 the same comparison with the ZEUS data [29] will be 



performed. In section 4 we shall finish with a summary and our conclusions. 



II. COMPARISON WITH RECENT HI DATA 



The recent HI data for diffractive photoproduction of dijets 30|, |3l| have several ad- 



vantages as compared to the earlier HI 



28l | and ZEUS [29] analyses. First, the integrated 



luminosity is three times higher than in the previous HI analysis [28] comparable to the 



luminosity in the ZEUS analysis [29]. Second, HI took data with \ow-E^ 



<:t 



30 



high-E^ 6 * 31 1 cuts, which allows the comparison with 28 1 and [29] . The exact two kine- 



31] and 
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TAB 
tion 



30 



: Kinematic cuts applied in the most recent HI analyses of dim-active dijet photoproduc- 



ail- 



HI low-E^ 6 * cuts 


HI high-f^f* cuts 


Q 2 < 0.01 GeV 2 


Q 2 < 0.01 GeV 2 


0.3 < y < 0.65 


0.3 < y < 0.65 


E ] T etl > 5 GeV 


M? w > 7.5 GeV 


E j T et2 > 4 GeV 


El et2 > 6.5 GeV 

_Z 


-1 < jfet\{2) < 2 


-1.5 < rf etl V) < 1.5 


zjp < 0.8 


zp < 1 


xjp < 0.03 


j;jp < 0.025 


\t\ < 1 GeV 2 


|t| < 1 GeV 2 


My < 1.6 GeV 


M Y < 1.6 GeV 



matic ranges are given in Tab. 1. These ranges for the low-E^ e * cuts are as in the previous 
HI analysis {^J and for the high-E^ 6 ' cuts are chosen as in the ZEUS analysis with two 
exceptions. In the ZEUS analysis the maximal cut on Q 2 is larger and the data are taken 
in an extended y range. The definition of the various variables can be found in the HI and 
ZEUS publications 28|, |29] and in our previous work 18|, [19] . Very important is the cut on 
xp. It is kept small in both analysis in order for the pomeron exchange to be dominant. 
We base our analysis on the \ow-Et data published in 30J, which differ from the data in 



3l| in the cut on zp influencing not only the experimental data, but also the NLO results 
in all variables except the distribution in zp. The preliminary low-E^ data 31 1 have been 



compared previously to our theoretical results at NLO in |32|]. In the experimental analysis 
as well as in the NLO calculatio ns, J ets are defined with the inclusive /cT-cluster algorithm 



with a distance parameter d — 1 



331 ] in the laboratory frame. At least two jets are required 



with the respective cuts on E^ 11 and E^ t2 , where E^ refers to the jet with the largest 
(second largest) E^ 1 . As is well known, the lower limits on the jet Et are chosen asymmetric 
in order to avoid an infrared sensitivity in those NLO cross section computations, which are 
integrated over E^* 34 1. 

Before we confront the calculated cross sections with the experimental data, we correct 
them for hadronization effects. The hadronization corrections are calculated by means of the 
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LO RAPGAP Monte Carlo generator 15J. The factors for the transformation of jets made 



up of stable hadrons to parton jets were supplied by the HI collaboration |30j, |3l|. Most of 
our calculations are done with the 'HI 2006 fit B' ^ DPDFs, since they give the smaller 
diffractive dijet cross sections as compared to the 'HI 2006 fit A'. The HI and ZEUS collab- 
orations constructed two more sets of DPDFs, which are called 'HI 2007 fit jets' and 'ZEUS 
DPDF SJ'. These fits are obtained through simultaneous fits to the diffractive inclusive and 
DIS dijet cross sections [35]. In these fits it is assumed that there is no factorization breaking 
in the diffractive DIS dijet cross sections. Including these cross sections in the fits leads to 
additional constraints, mostly for the diffractive gluon distribution. On average the 'HI 2007 
fit jets' is similar to the 'HI 2006 fit B' except for the gluon distribution at large momentum 
fraction and small factorization scale. In the following analysis we shall disregard these new 
DPDF sets, since they would be compatible with the factorization test of the photoproduc- 
tion data only, if we restricted these tests to the case that the resolved part has the breaking 
and not the direct part, which has the same theoretical structure as the DIS dijet cross 



section. Results with the 'HI 2007 fit jets' can be found in [3QJ, |3jJ. The 'HI 2006 fits A,B' 
are based on n/ = 3 massless flavors. The production of charm quarks was treated in the 
Fixed-Flavor Number Scheme (FFNS) in NLO with non-zero charm-quark mass yielding a 
diffractive F$- This _F 2 C is contained in the 'HI 2006 fits A,B' parameterizations and is then 
converted by us into a charm PDF using the LO expression F% (x, Q 2 ) = 2xe 2 f c (x, Q 2 ), where 
e c = 2/3 is the electric charge of the charm quark. The bottom contribution was neglected. 
This assumption simplifies the calculations considerably. Since the charm contribution from 
the Pomeron is small, this should be a good approximation. We then take rif — 4 with 
A|| = 0.347 GeV, which corresponds to the value used in the DPDFs 'HI 2006 fits A,B' 
[4j. For the pho ton PDF we have chosen the NLO GRV parametrization transformed to the 



MS scheme 



36]. 



As it is clear from the discussion of the various preliminary analyses of the HI and ZEUS 
collaborations, there are two questions which we would like to answer from the comparison 
with the recent HI and the ZEUS data. The first question is whether a suppression factor 
(sometimes also called rapidity-gap survival probability), which differs substantially from 
one, is needed to describe the data. The second question is whether the data are also 
consistent with a suppression factor applied to the resolved cross section only. To give an 
answer to these two questions we calculated first the cross sections with no suppression 
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factor (R = 1 in the following figures) with a theoretical error obtained from varying the 
common scale of renormalization and factorization by factors of 0.5 and 2 around the default 
value (highest E 3 ^ 1 ). In a second step we show the results for the same differential cross 
sections with a global suppression factor, adjusted to da/dE 3 ^ 11 at the smallest E^ a -hm. 
As in the experimental analyses 30|, 131] , we consider the differential cross sections in the 
variables x° fes , zp, log 10 (xj>)> E 3 T etl , M x , M 12 ^% \Ar]^ ts \ and W. The definitionof the 
variables is given in the experimental papers 



28H3lj or in our earlier work 



18 



,1191,1371. In 



the latter references also the relevant formulas for the calculation of the dijet cross sections 
can be found. 

For the low-E^ cuts, the resulting suppression factor is R = 0.50 ± 0.09, which gives in 
the lowest E^ n -hm. a cross section equal to the experimental data point. The error comes 
from the combined experimental statistical and systematic error. The theoretical error due 
to the scale variation is taken into account when comparing to the various distributions. 
The result of this comparison is shown in Figs. 3a-i. With the exception of Figs. 3d and 3f, 
where the comparisons of da/dE 3 ^ 1 and dajdM\2 are shown, all other plots are such that 
the data points lie outside the error band based on the scale variation for the unsuppressed 
case. However, the predictions with suppression R = 0.50 agree nicely with the data inside 
the error bands from the scale variation. Most of the data points even agree with the 
R = 0.50 predictions inside the much smaller experimental errors. The M^-distribution 
agrees only for the second bin. The reason for the disagreement of the two other Mx-bins 
might be that in the theoretical results the variable Mx is defined without the remnant 
contributions, which, however, are taken into account in the experimental definition of Mx- 
Further exceptions are the cross sections da/dE 3 ^ 11 and da/dMi 2 , which are related. In 
da/dE 3 ^ 11 (see Fig. 3d) the predictions for the second and third bins lie practically outside 
the data points with their errors (note the logarithmic scale). For R = 1 and R = 0.50 the 
cross sections fall off stronger with increasing E 3 ^ 1 than the data, the normalization being 
of course about two times larger for R — 1. In particular, the third data point agrees almost 
with the R = 1 prediction. This means that the suppression decreases with increasing Ep 11 . 
This behavior was already apparent when we analyzed the first preliminary HI data 18l. Il9|. 
Such a behavior points in the direction that a suppression of the resolved cross section only 
would give better agreement with the data, as we shall see below. The same observations can 
be made by looking at da/dMw in Fig. 3f. The survival probability R = 0.50 ± 0.09 agrees 
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FIG. 3: Differential cross sections for diffractive dijet photoproduction as measured by HI with 
low-E^ cuts and compared to NLO QCD without (R = 1) and with (R = 0.50) global suppression 
(color online). 



with the result in [30(], which quotes R = 0.58 ± 0.01 (stat.) ± 0.12 (syst.), determined 
by fitting the integrated cross section. From our comparison we conclude that the low-E^ 
data show a global suppression of the order of two in complete agreement with the results 



in 



1<< 



19| and [28| based on earlier preliminary [l^] and final HI data 28] 



Next we want to answer the second question, whether the data could be consistent with 
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a suppression of the resolved component only, whose definition is not unique, but rather 
factorization scale and scheme dependent. For this purpose we have calculated the cross 
sections in two additional versions: (i) suppression of the resolved cross section in the MS 
scheme and (ii) suppression of this resolved cross section plus that part of the NLO direct part 



which depends on the factorization scale at the photon vertex c 



Of course, the needed 



suppression factors for the two versions will be different. We determine the suppression 
factors again by fitting the measured da/dE 3 ^ 11 for the lowest E^-hm (see Fig. 4d). 
Then, the suppression factor for version (i) is R = 0.40 (denoted res in the figures), and 
for version (ii) it is R = 0.37 (denoted res+dir-IS). The comparison with the HI data of 
da/dx° bs , da/dz$ s , da/dlog 10 (x P ), da/dE ] T etl , da/dM x , da/dM 12 , da/d^ ets , da/d\A^ ets \ 
and da/dW is shown in Figs. 4a-i, where we have also plotted the prediction for the global 
suppression (direct and resolved) with R = 0.50, already shown in Figs. 3a-i. Looking at the 
Figs. 4a-i we can distinguish three groups of results from the comparison with the data. In 
the first group, the cross sections as functions of z^ s , /ogio(xp), M i2 , ff ets , lAr/p 6 * 5 and W, 
the agreement with the global suppression (R = 0.50) and the resolved suppression (R = 0.40 
or R = 0.37) is comparable. In the second group, namely for da/dE^ n , the agreement is 
better for the resolved suppression only. In the third group, da/dx° bs and da/dMx, the 
agreement with the resolved suppression is worse than with the global suppression. In 
particular, for da/dx° bs , which is usually considered as the characteristic distribution for 
distinguishing global versus resolved suppression, the agreement with resolved suppression 
does not improve. Unfortunately, this cross section has the largest hadronic corrections 
of the order of (20 — 30)% j^]. Second, also for the usual photoproduction of dijets the 
comparison between data and theoretical results has similar problems in the large x° bs -bin 
39], although the E^ 1 cut is much larger there. On the other hand, for the cross sections 



da/dE^ n (and da/dMi 2 ) the agreement improves considerably (and somewhat) with the 
suppression of the resolved part only (note the logarithmic scale in Fig. 4d). Here, of course, 
we must admit that the suppression factor could be Et- dependent, although we see no 
theoretical reason for such a dependence. 

We also checked for two distributions whether the predictions for resolved suppression 
depend on the chosen diffractive PDFs. For this purpose we have calculated for the two 
cases da/dzp and da/dE J T etl the cross sections with the 'HI 2006 fit A' parton distributions 
4j. The results are compared in Figs. 5a and b to the results with the 'HI 2006 fit B' and the 
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FIG. 4: Differential cross sections for diffractive dijet photoproduction as measured by HI with low- 
Ej? 1 cuts and compared to NLO QCD with global, resolved, and resolved/direct-IS suppression. 
Note that some of the theoretical predictions coincide with the experimental values. 

experimental data. Of course, since the 'HI 2006 fit A' PDFs have a larger gluon component 
at large z, the cross sections are larger and therefore need a larger suppression of R = 0.32. 
Note that in the published low-i^ e * HI analysis as well as in the comparison presented here 
the contribution from the largest has been removed from all other distributions. 

From Figs. 5a and b we conclude that the dependence on the chosen DPDFs is then weaker, 



15 



1000 



800 



•8 600 



■§ 400 



200 



ep — > e+2jets+X'+Y 



ep — > e+2jets+X'+Y 



-J — 1 HI 99-00 

NLO, R=0.40 (res) (HI 2006 Fit B) 
NLO, R=0.32 (res) (HI 2006 Fit A) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

obs 



10' 



o 

r 

Q. 



w 

I 10 r > 

T3 



-S — i HI 99-00 

NLO, R=0.40 (res) (HI 2006 Fit B) 
NLO, R=0.32 (res) (HI 2006 Fit A) 



9 10 11 12 13 14 

jetl 

Ep [GeV] 



FIG. 5: Differential cross sections for diffractive dijet photoproduction as measured by HI with 



\o~w-Ej? 1 cuts and compared to NLO QCD with resolved suppression and two different DPDFs. 



but that 'HI 2006 fit B' is still favored over 'HI 2006 fit A'. In total, we are tempted to 
conclude from the comparisons in Figs. 4a-i that the predictions with a resolved-only (or 
resolved+direct-IS) suppression are consistent with the new low-£f* HI data |30|. 

The same comparison of the high-E^ 6 ' data of HI [3l| with the various theoretical pre- 
dictions is shown in the following figures. The global suppression factor is obtained again 
from a fit to the smallest E^ a -hm.. It is equal to R = 0.62 ± 0.16, again in agreement with 
the HI result R = 0.62 ± 0.03 (stat.) ± 0.14 (syst.) [3l| obtained with our theoretical 
cross sections. The comparisons of the same cross sections as in the low-E^ 1 comparison 
are shown in Figs. 6a-i for the two cases R = 1 (no suppression) and R = 0.62 (global sup- 
pression). As before with the exception of da/dE^ 1 and dojdMyi in Fig. 6d and Fig. 6f, 
most of the data points lie outside the R = 1 results with their error bands and agree with 
the suppressed prediction with R = 0.62 inside the respective errors. However, compared 
to the results in Figs. 3a-i the distinction between the R = 1 band and the R = 0.62 band 
and the data is somewhat less pronounced. We also tested the prediction for the resolved 
(resolved+direct-IS) suppression, which is shown in Figs. 7a-i. The suppression factor fitted 
to the smallest E^ a -bm came out as R = 0.38 (res) and R = 0.30 (res+dir-IS). In most 
of the comparisons it is hard to observe any preference for the direct plus resolved (global) 
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FIG. 6: Differential cross sections for diffractive dijet photoproduction as measured by HI with 
high-E^ 6 * cuts and compared to NLO QCD without (R = 1) and with (R = 0.62) global suppression 
(color online). 

suppression against the resolved suppression only. We remark that the suppression factor 
for the global suppression is increased by 24%, if we go from the low-£^ et to the high-i?jf* 
data, whereas for the resolved suppression the difference is only 5%. Under the assumption 
that the suppression factor should not depend on E^ 11 , we would conclude that the version 
with the resolved suppression would be preferred. In Figs. 8a and b we tested the resolved 
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FIG. 7: Differential cross sections for diffractive dijet photoproduction as measured by HI with 
high-.Ejf' cuts and compared to NLO QCD with global, resolved, and resolved/direct-IS suppres- 
sion. 

suppression model against the choice of the two DPDFs, fit A versus fit B, with the result 
that this dependence is weak if we adjust the suppression factor, which is R = 0.16 for 
the 'HI 2006 fit A'. The general conclusions from the high-E^ 6 * comparison are very much 
the same as from the low-E^ 6 * comparison. A global suppression is definitely observed and 
the version with resolved suppression explains the data almost as well as with the global 
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FIG. 8: Differential cross sections for diffractive dijet photoproduction as measured by HI with 
high-i?;jf* cuts and compared to NLO QCD with resolved suppression and two different DPDFs. 



suppression. 



III. COMPARISON WITH ZEUS DATA 



In this section we shall compare our 



which was published just recently 



predictions with the final analysis of the ZEUS data, 



291 ]. in order to see whether they are consistent with 
the large- E 3 ^ data of HI. The kinematic cuts are almost the same as in the high-£^ et HI 
measurements. They are given in Tab. 2. The only major difference to the HI cuts in Tab. 
1 is the larger range in the variable y. Therefore the ZEUS cross sections will be larger than 
the corresponding HI cross sections. The different cuts on Q 2 and \t\ have little influence. 
For example, the larger |t|-cut in Tab. 2 as compared to Tab. 1 increases the cross section 
only by 0.2%. The constraint on My is not explicitly given in the ZEUS publication 29 ]. 
They give the cross section for the case that the diffractive final state consists only of the 
proton. For this they correct their measured cross section by subtracting in all bins the 
estimated contribution of a proton-dissociative background of 16%. When comparing to the 
theoretical predictions they do the reverse and multiply the cross section with the factor 
0.87, in order to correct for the proton-dissociative contributions, which are contained in the 
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TABLE 
duction 



I: Kinematic cuts applied in the most recent ZEUS analysis of diffractive dijet photopro- 



29]. 



ZEUS cuts 



Q 2 < 1 GeV 2 
0.2 < y < 0.85 
E J T etl > 7.5 GeV 
Eg t2 > 6.5 GeV 
-1.5 < rf etl W < 1.5 
xjp < 0.025 
\t\ < 5 GeV 2 



DPDFs 'HI 2006 fit A' and 'HI 2006 fit B' by requiring M Y < 1.6 GeV. We do not follow 
this procedure. Instead we leave the theoretical cross sections unchanged, i.e. they contain a 
proton-dissociative contribution with My < 1.6 GeV, and multiply the ZEUS cross sections 
by 1.15 to include the proton-dissociative contribution. Since the ZEUS collaboration did 
measurements only for the high-Z^ 6 * cuts, E 3 ^ > 7.5 (6.5) GeV, we can only compare 
to those. In this comparison we shall follow the same strategy as before. We first compare 
to the predictions with no suppression (R = 1) and then determine a suppression factor 
by fitting da/dE^ 1 to the smallest Ej? -bin. Then we compare to the cross sections as a 
function of the seven observables x° bs , Zp s , xjp, E^ 1 , y, Mx and "if etx instead of the nine 
variables in the HI analysis. The distribution in y is equivalent to the H^-distribution in 



311 ] . The theoretical predictions for these differential cross sections with no suppression 
factor (R = 1) are shown in Figs. 9a-g, together with their scale errors and compared to 
the ZEUS data points. Except for the x° bs - and .^^-distributions, most of the data points 
lie outside the theoretical error bands for R — 1. In particular, in Figs. 9b, c, e, f and g, 
2, 3, 4, 4 and 5 points lie outside. This means that most of the data points disagree with 
the unsuppressed prediction. Next, we determine the suppression factor from the measured 
da/dE 3 T etl at the lowest E J T etl -bin, 7.5 GeV < E J T etl < 9.5 GeV, and obtain R = 0.71. As 
a curiosity, we remark that this factor is larger by a factor of 1.15 than the suppression 
factor from the analysis of the high-i?^ e< data from HI. This factor is exactly equal to the 
correction factor we had to apply to restore the dissociative proton contribution. Without 
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FIG. 9: Differential cross sections for diffractive dijet photoproduction as measured by ZEUS and 
compared to NLO QCD without (R = 1) and with (R = 0.71) global suppression (color online). 

this correction factor the suppression factor following from the ZEUS analysis would be in 
perfect agreement with the factor in the HI analysis. Taking the total experimental error of 
±7% from the experimental cross section da/dE^ 11 in the first bin into account, the ZEUS 
suppression factor is 0.71 ± 0.05 to be compared to 0.62 ± 0.14 in the HI analysis [3 If . so 
that both suppression factors agree inside the experimental errors. 

If we now check how the predictions for R = 0.71 compare to the data points inside 
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the theoretical errors, we observe from Figs. 9a-g that with the exception of da/dzp 8 and 
da/dE^ n the majority of the data points agree with the predictions. This is quite con- 
sistent with the HI analysis, discussed in the previous section, and leads to the conclusion 
that also the ZEUS data agree much better with the suppressed predictions than with the 
unsuppressed prediction. In particular, the global suppression factor agrees with the global 
suppression factor obtained from the analysis of the HI data inside the experimental error. 

Similarly as in the previous section we compared the ZEUS data also with the assumption 
that the suppression results only from the resolved cross section. Here we consider again the 
two versions: (i) only resolved suppression (res) and (ii) resolved plus direct suppression of 
the initial-state singular part (res+dir-IS). For these two models we obtain the suppression 
factors R = 0.53 and R = 0.45, respectively, where these suppression factors are again 
obtained by fitting the data point at the first bin of da/dE 3 ^ 11 . The comparison to the 
global suppression with R = 0.71 and to the data is shown in Figs. lOa-g. In general, we 
observe that the difference between global suppression and resolved suppression is small, i.e. 
the data points agree with the resolved suppression as well as with the global suppression. 

In Figs. 11a and b the difference between 'HI 2006 fit B' and 'HI 2006 fit A' is shown 
again for the case of the resolved suppression. In both figures we observe that the fit A 
suppression with the suppression factor R = 0.27 agrees better with the data than with the 
factor R = 0.53 for the fit B suppression. In particular, for da/dE^ n the agreement with 
the three data points is perfect (note the logarithmic scale). 

In our analysis of the ZEUS data so far we assumed that the measurements with the large 
rapidity gap (LRG) method of ZEUS in 29| are such that with this method the same inclusive 
diffractive DIS cross section is measured as in in the HI measurement of this cross section 
with the LRG method [4], on which the fits of the DPDFs 'HI 2006 fits A,B' are based. This 
is actually not true, and this problem has been analysed by the ZEUS collaboration in their 
publication, in which they present their data for the inclusive diffractive DIS cross sections 
using different definitions for these cross sections 40j . They find that their LRG cross section 
has to be corrected by two factors in order to make it agree with the HI measurement of 
the diffractive DIS cross section and with the prediction of this cross section based on the 
'HI 2006 fit B'. First, a factor of 0.91 ± 0.07 was estimated with the PYTHIA Monte Carlo, 
so that the ZEUS cross sections correspond to a dissociative background with My < 1.6 
GeV. Second, even then the so corrected ZEUS diffractive DIS cross section was still 13% 
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FIG. 10: Differential cross sections for diffractive dijet photoproduction as measured by ZEUS and 
compared to NLO QCD with global, resolved, and resolved/direct-IS suppression. 



larger than the corresponding HI diffractive DIS cross section. This amounts to a total 
correction factor of 0.79 ± 0.06. If we apply this correction factor to our ZEUS suppression 
factors, we obtain a global (resolved-only) suppression of 0.56 ± 0.05 rather than 0.71 (0.42 
± 0.04 rather than 0.53), i.e. suppression factors which are closer to the results found for 
the similar high-i^* HI analysis. Here, the errors refer only to the errors coming from 
the renormalization of the ZEUS data and do not include the experimental stat./syst. and 
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FIG. 11: Differential cross sections for diffractive dijet photoproduction as measured by ZEUS and 
compared to NLO QCD with resolved suppression and two different DPDFs. 



theoretical scale errors. 



IV. CONCLUSION 



In this paper, we confronted the final HERA data on the diffractive photoproduction of 
dijets, as published by the ZEUS [29] and HI [30, |3l| collaborations, with our NLO QCD 
calculations in order to see if factorization breaking effects in the resolved and eventually 
direct cross sections can be established consistently from both data sets. The new comparison 
is even more conclusive than the one published previously in our invited review as the 
proton beam energy of 920 GeV, the cuts on the jet transverse energies oiE J T etl{2) > 7.5 (6.5) 
GeV, and the cut on the momentum fraction carried by the diffractive exchange xp < 0.025 
are now the same in both experiments, which was not the case before. At the same time, 
some experimental cuts are still different. In particular, the momentum fraction y transferred 
by the electron to the hadronic system is larger for ZEUS than for HI. We also re-computed 
in NLO QCD the cross sections for the lower cuts on E^ 11 ^ > 5 (4) GeV, which have 



been re-measured by the HI collaboration 
their previous low-luminosity data set 



30 



31| in order to establish consistency with the 



We found that the large majority of HI and 
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TABLE III: Suppression factors for global and resolved-only suppression in the low-E^ and high- 
E^* analyses of HI and the ZEUS analysis before and after rescaling their data by a factor of 
1.15. 



Suppression factor 


HI low-Ef * 


HI high-^ et 


ZEUS 


ZEUS renormalized 


global 


0.50 


0.62 


0.71 


0.56 ± 0.05 


resolved-only 


0.40 


0.38 


0.53 


0.42 ± 0.04 


res+dir-IS 


0.37 


0.30 


0.45 


0.36 ± 0.03 


res, HI 2006 fit A 


0.32 


0.16 


0.27 


0.21 ± 0.01 



ZEUS data points lay below the NLO QCD predictions, even when using the 'HI 2006 fit 
B' diffractive PDFs with small gluon density at large fractional momentum and taking into 
account the experimental (statistical and systematic) and theoretical (scale variation) errors. 



The data at larger E^ (or M 12 ) tended to agree better with the NLO QCD predictions 
than those at small E^ a . 

By fitting the lowest (and dominant) bins in the three -E^ efl -distributions, we established 
the amounts by which both the direct and resolved NLO QCD cross sections had to be 
reduced to find agreement with the data. These suppression factors are shown in the second 
line ('global') of Tab. 3 for the low-E^ 6 * and high-E^ 6 * analyses of HI, the ZEUS analysis 
(multiplied by a factor of 1.15 to allow for proton dissociation), and the ZEUS analysis 
renormalized by a factor of 0.79 ± 0.06 for My < 1.6 GeV and correspondance with the 
HI measurements and DPDF fits [40]. The first, second and third factors agreed very well 
with those found by the experimental collaborations when fitting multiple distributions or 
total cross sections. The fourth factor ('ZEUS renormalized') agrees better with the second 
factor, relevant for the similar high-i^ 6 * HI analysis. We also tested the hypotheses that 
factorization breaking is only present in the resolved (third line) or the resolved and the 
related initial-state singular part of the direct photoproduction cross sections (fourth line). 
Both hypotheses gave very similar results and described the data sets almost as well as the 
predictions with global factorization breaking. The suppression factors applicable to just 
;he resolved cross section are shown in the third line of Tab. 3. As observed previously 
19] ]. they agree very well with absorptive-model predictions jl3 |. We conclude that in this 
case the suppression factors do not show a significant E? -dependence, in particular when 
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renormalizing the ZEUS data as described above. The fact that no E^ 1 - dependence is 
visible here can, of course, be explained by the fact that the resolved cross section falls more 
steeply with E^ 1 than the direct one 12|. Finally, we investigated whether these conclusions 
depended on the diffractive PDFs by comparing the results with resolved-only suppression of 
the 'HI 2006 fit B' to those obtained with the 'HI 2006 fit A' (last line in Tab. 3). Since the 
latter has a larger gluon density at large momentum fraction, the suppression had to be more 
important. The fit A results then tended to describe the high-.Ejf* HI data and the ZEUS 
data slightly better, in particular in the ^^-distribution, which should be directly sensitive 
to the DPDFs, but the low-£^ et HI data slightly worse. Unfortunately the experimental 
and theoretical errors are still too large to draw any strong conclusions. 

While the epoch of HERA experiments has now ended and an International Linear Col- 
lider may not be built in the near future, it will be very interesting to investigate diffractive 
physics at the LHC Suprisingly, proton-proton and heavy-ion collisions at the LHC can also 
be a source of high-energy photon collisions, and this may open up a whole new field of 
investigation for diffractive dijet photoproduction |41 |. 
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